Spontaneous inhomogeneous phases in ultracold dipolar Fermi gases 
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We study the collapse of ultracold fermionic gases into inhomogeneous states due to strong dipolar 
interaction in both 2D and 3D. Depending on the dimensionality, we find that two different types 
of inhomogeneous states are stabilized once the dipole moment reaches a critical value d > d c : 
the stripe phase and phase separation between high and low densities. In 2D, we prove that the 
stripe phase is always favored for d > d c , regardless of the microscopic details of the system. In 
3D, the one-loop perturbative calculation suggests that the same type of instability leads to phase 
separation. Experimental detection and finite-temperature effects are discussed. 
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I. INTRODUCTION 

The rapid experimental development of ultracold 
atomic and molecular physics has opened up new op- 
portunities to study quantum many-body systems with 
electric and magnetic dipolar interactions [lH9(. An im- 
portant feature of dipolar interaction is its explicit spa- 
tial anisotropy of the d,.2_ 3z 2-type when dipole moments 
are aligned by external fields. For the fermionic dipo- 
lar systems, 40 K- 87 Rb has been cooled down almost to 
quantum-degeneracy Anisotropic Fermi liquid theo- 
ries of the single particle and collective properties have 
been investigated 

EMI- 

Furthermore, unconventional 
Cooper pairing structures have been studied, including 
the TV-channel pairing in the single component systems 
p7H22| , and the competition between the s + d-wave sin- 
glet and the p z -wave triplet channels [23T - I25T ] . In par- 
ticular, a novel pairing state of the s + ip t ype broken 
time-reversal symmetry has been pointed out |25| . More- 
over, magnetic dipolar systems of fermions have also been 
experimentally realized @. Exotic states of the ferro- 
nematic Fermi liquid and unconventional magnetic states 
have been predicted [H, H3] ■ 

Fermionic systems can spontaneously break transla- 
tional symmetry in both charge and spin channels. Many 
years ago, Overhauser pointed out that even in the weak 
coupling regime of the interacting electron gas, a 2k f 
spin-density wave state always wins over the uniform 
paramagnetic state at the Hartree-Fock level [28j . How- 
ever, correlation effects may suppress this instability and 
such a state has not been experimentally confirmed. Re- 
cently, quantum liquid crystal phases in strongly corre- 
lated systems have been intensively studied, particularly 
in doped Mott insulators [29[. Stripe ordering has been 
observed in high T c cuprates, other transition metal ox- 
ides, and quantum Hall systems at high Landau levels 

[1113. 

In this article, we study the instability toward the 
spontaneous inhomogeneous phase in dipolar fermionic 
systems. In two dimensional (2D) systems, the strongest 
density-channel instability occurs at non-zero momen- 



tum, which drives the density wave states under strong 
dipolar interactions. This effect is based on the pecu- 
liar feature of the Fourier transform of the dipolar in- 
teraction, which is robust against microscopic details. 
However in three dimensions (3D), the instability at the 
one-loop level occurs at zero momentum, thus it leads to 
phase separation into high and low density regions. 

This paper is organized as follows: In Sec. HH we 
construct the model Hamiltonian for a single-component 
dipolar Fermi gas and present the calculation of the static 
susceptibility of density fluctuations via diagrammatic 
expansion. Then, we study the instability of the homo- 
geneous ground state in 2D and 3D in Sees. HH and Hvl 
correspondingly. Here, we present both a nonperturba- 
tive conclusion, which is based on the analytic property 
of the susceptibility and the leading-order perturbative 
calculation, which is use to testify the nonperturbativc 
conclusion. The experimental detection is then discussed 
in Sec. [V] Finally, we conclude our paper by studying 
the contributions of thermal fluctuations at finite tem- 
perature and the effects of higher-order terms in Sec. I VII 



II. MODEL 

We consider the single-component dipolar Fermi gas 
with dipole moment aligned by an external electric field. 
The long-distance physics of this system is described by 
the following Hamiltonian 

k k,k f ,q 

with eg being the energy dispersion relation and /i the 
chemical potential. In 2D, the dipolar interaction in the 
momentum space V2d(q) takes the form (l3| . 

V 2 D{q) = 2TTd 2 P 2 {cos0)[- -q] + irqd 2 sm 2 9 cos 2(f> q , (2) 

where d is the dipole moment; 4> q is the azimuthal angle 
of the momentum q and 9 is the angle between the di- 
rection of the dipoles and the z-axis. e is a short-range 
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FIG. 1: (online color) Diagrammatic expansion for Ti^'(q). 
The solid lines are free-fermion propagators and the dashed 
lines represent the dipolar interaction. 




cutoff roughly equal to the thickness of the system along 
the z-direction. Equation ([2]) is valid for qe <C 1. How- 
ever, as shown below, the main conclusion of this paper 
only relies on the value of V(q) at small q, and hence 
all the qualitative results remain invariant even if more 
accurate and more complicated interactions were taken 
into account. We shall also emphasize here that for 6 = 0, 
V2d(<z) is invariant under SO (2) rotations. However, for 
any 0^0, the rotational symmetry is reduced to two-fold 
(C2). This symmetry property is crucial as a system un- 
dergoes the transition to a charge-density-wave (CDW) 
phase. As will be discussed below, in a system with SO(2) 
symmetry, the CDW instability usually leads to the for- 
mation of a triangular lattice [33| (a 2D Wigner crystal). 
On the other hand, the same instability in a system with 
C2 symmetry is compatible with the formation of a uni- 
directional CDW, i.e., a stripe phase [34l |. 

Notice that Eq. @ contains an isotropic part which 
depends on the microscopic cut-off e and an anisotropic 
part which depends on the azimuthal angle of <f> q . We 
first consider the purely anisotropic case at 

9 = 9 Q = cos- 1 -^, (3) 

and then the interaction simplifies into 

ViT>{q) = ^pzcos2</v (4) 

The major feature of Eq. ((4]) is its linear dependence on 
q, which play important roles in driving the instability at 
non-zero wavevectors. 

The stability of a system against density fluctuations 
is determined by the static susceptibility of density fluc- 
tuations: 

n 2D («7> = 0) = (p(q,cj = 0)p(-q,uj = 0)), (5) 

where 

fc,0 

is the density operator. The inverse of I12d(<T, w = 0) 
is the energy gap A 2 d(<z) of creating a density wave at 
momentum q. Using the technique of the diagrammatic 
expansion, A2d(<?) can be computed order by order as 

A 2D (q) = n 2D (q> = 0)- 1 

= in™(<?,u; = o)}- 1 + vMq), (?) 



FIG. 2: (online color) The stripe phase in 2D dipolar gases. 
The three dotted arrows mark the x, y and z axes. The lighter 
(darker) regions in the 2D plane represent higher (lower) den- 
sity. The solid line with an arrow represents the direction of 
the external field, which aligns the direction of the dipoles. 



where H im (q,ui) is the one-particle- irreducible density- 
density correlation [3^, H|| shown in Fig. [TJ The sta- 
bility criteria of the homogeneous phase corresponds to 
the energy gap A2d(<7) being positive-definite for any q. 
Otherwise, the creation of a density wave could lower the 
total energy of the system leading to the condensation of 
density waves. Obviously, A2d(<7) is positive-define for 
any q in the noninteracting Fermi gas. This stability re- 
mains for a weakly-interacting Fermi liquid. On the other 
hand, if the interaction has an attractive channel [e.g. the 
dipolar interaction in Eq. an inhomogeneous state 

could take advantage of the attractive interaction. Under 
strong enough attractive interactions, i.e., V2D —00, 
the homogeneous ground state is unstable. 



III. STRIPE PHASES IN 2D 

We begin from the weak coupling side, by tuning up 
the interaction strength. In general, there is a critical 
dipole moment d c such that A2d(<?) > at any q for 
d < d c , but A 2 d(<?) < at some q for d > d c . At 
d = d c , A2d(±Q) = for the momenta ±Q and remain 
positive at other momenta [37| . If we increase d further 
above d c , the density wave fluctuations with momentum 
q ~ Q undergo an instability and condense. Depending 
on whether Q is or not, this instability has two differ- 
ent fates. If Q > 0, the condensation carries a nonzero 
momentum, which leads to a stripe state: a unidirec- 
tional charge-density-wave state with wavevector Q. On 
the other hand, if Q = 0, the condensation takes place 
at an infinite wavelength and results in phase separation 
between high and low density regions. Here, the typical 
size of these high (low) density regions and their spatial 
arrangements are determined by the microscopic details 
of the system. In the particular case of a dipolar gas, the 
phase separation is sensitive to the details of the short- 
range behaviors of the interaction. 
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Due to the 2D space-inversion symmetry r — > —r, 
PaoHSi w )] 1 i s an even function of q. As a result, the 
Taylor expansion of A 2 d(<z) at small momentum must 
take the following form 



A 2D (<7) =[^n(Q = 0) 



2nd 2 



-q cos 2(j) q 



(8) 



At small momentum, q ~ 0, Eq. ([5]) was dominated by 
the first two terms. And hence the point q = is a local 
saddle point of A 2 d(<z') (i.e., a local minimum along the 
line of 4> q = and a local maximum along <f> q = tt/2). 
Since g = is a saddle point, instead of the minimum 
point of A 2 d (<?*), as d increases, A 2 d(q) at some finite mo- 
mentum will become negative before A 2 D(<f = 0) could 
reach zero. Thus, Q must be a nonzero value, which in- 
dicates that the system will collapse into a stripe phase 
for d > d c , instead of phase separating. This is the main 
conclusion of this paper. It is noteworthy that this con- 
clusion only depends on the small-momentum (i.e. long- 
range) behaviors of the interaction V(q). Therefore, the 
stripe instability is a universal property of the 2D dipolar 
Fermi gas at d > d c , insensitive to the microscopic details 
of the system. In 2D electron g similar conclusion 

has been found in the classical limit at high temperature 
when the quantum fluctuations are negligible 38] . 

We emphasize that the conclusion of Q > is a non- 
perturbative result which remains valid to any order in 
the loop expansion of Fig. [TJ The only assumption re- 
quired here is the analyticity of ni^if, oj — 0) as a func- 
tion of q near q = 0. In a Fermi liquid, it is well-known 
that H im (q,u) — 0) is non-analytic at q = 2fcp where fcp 
is the Fermi wavevector. However, there is no reason of a 
non-analytic behavior at small momentum in any Fermi 
liquid to our knowledge. The higher-order diagrams in 
Fig. [T] contain interaction lines which are non-analytic 
at small momenta. However, due to the fact that dipo- 
lar interactions are short-ranged in 2D, we do not expect 
any singularity for n2™(<f, oj — 0) around q = 0. 

We check this conclusion using the loop expansion in 
Fig. Q] to the one-loop level for a 2D system at T = 
in continuum (e^ = k 2 /2m with m being the mass of the 
particles) and find 




FIG. 3: (online color) One-loop calculation of A(q) in (a) 2D 
and (b) 3D. In the 2D case, we choose 4>q = § an d for the 
3D, 8q = 7r/2. From top to bottom, the dipole momentum for 
each of the curves are: d/d c = 0, 0.9, 1, and 1.1 respectively. 
The dashed vertical line marks the momentum 2/cf- 



d. There, the density need to be higher than a criti- 
cal value to reach the stripe phase. For KRb (d — 0.57 
D) the stripe phase requires the interparticle spacing to 
be smaller than ~ 1CP 4 cm, and for LiCs (d ~ 5.5 D), 
the critical spacing is ~ 10~ 2 cm. Wc emphasize that 
Q and d c we found here are based on the leading order 
approximation, and can be renormalized by higher order 
corrections. However, the fact that Q > remains valid, 
which is one of the key conclusion of this work. 



For ^ 9q, the same conclusion applies for the single 
component case. The e-term in Eq. ([5]) is momentum 
independent, which corresponds to short-range contact 
interaction. Thus it has no effect on the single component 
fermions, and hence 



V(q) = qird 2 [-2P 2 (cos6) + sin 2 6 cos 20, 



(10) 



A 2D (<?) 



2 -\ 

771 



1 



1 



0^7(1 -a) 



dl 



a;cos20g-], (9) 



where x = q/2kF and d c is the critical dipole momentum 
d c = \Ji/ (2kpm). <fiq~ is the azimuthal angle of the mo- 
mentum q and rj(x) is a step function with r)(x) = 1 for 
x > and rj(x) = for x < 0. As shown in Fig. 3(a) 
the order wavevector Q — 2kp. In experiments, it is of- 
ten more convenient to vary the density for a constant 



Along cj> q = | 
strongest instability. 



which remains the direction of the 



V(q) = -2qird 2 cos 2 6, 



(11) 



which is negative and grows linearly with q for an arbi- 
trary value of 6. Thus, the above conclusion still holds. 
We notice that even for 8 = 0, i.e., when the dipoles are 
perpendicular to the plane, the stripe instability still ex- 
ists at large values of d although the interaction is now 



4 



repulsive. For such an isotropic system (at 6 = 0), it is 
often assumed that three charge density waves form a tri- 
angular lattice to minimize the breaking of the rotational 
symmetry [33]. However, if the rotational symmetry is 
spontaneously broken before the homogeneous state be- 
comes unstable, a stripe phase can be stabilized at zero 
temperature [341 ] . In this case, the stripe ordering spon- 
taneously breaks both the rotational symmetry and the 
translational symmetry in the direction perpendicular to 
the stripes, while only the latter is spontaneously break- 
ing in the stripe phases at 9 ^ 0. Due to its unique 
symmetry breaking pattern, the low-energy properties of 
the stripe phase at 8 — is described by the quantum 
McMillan-de Gennes theory, which is fundamentally dif- 
ferent from other similar phases at 9 ^ [34[ • As a result, 
the stripe phase at 9 = is unstable at any finite tem- 
perature due to thermal fluctuations [3~il l39j , while the 
stripe phases at 9 7^ has a power-law quasi-long range 
order blow a transition temperature Tkt (See below). 



where the spin indices a and a' run from —N/2 to +N/2 
for some integer N. Here, V{q) remains the same as 
the single-component case in Eqs. ([2]) and (fT2|) . For 
a single-component dipolar gas, due to Pauli exclusion 
principle, the e-term in Eq. ([2]) (short-range interactions) 
has no contribution to any physical properties of the sys- 
tem. However, for multi-component Fermi gases, this 
term naturally exists in the inter-component interactions 
and suppresses the stripe instability if it is repulsive. For 
simplicity we only consider the case of N = 2. For 9 > 9q 
(9 < 6 ), this short-range part is attractive (repulsive), 
and hence the stripe state become unfavorable for 9 < 9q. 
Furthermore, such an finite momentum instability does 
not exist in the spin channel because the spin channel 
interaction arises from the exchange interaction which is 
regular for small vector of q. 



V. EXPERIMENTAL DETECTION 



IV. FERMIONS IN 3D AND 
MULTICOMPONENT FERMIONS 

In this section, we study two more classes of systems 
with similar dipolar interactions: (a) dipolar gases in 3D 
and (b) 2D dipolar gases composed by multicomponent 
fermions. 

In 3D, the dipolar interaction is 

V 3 B(q) = ^fp2(cos9 q ~), (12) 

with 9 q measuring the angle between q and the dipole 
moment. Due to the structure of V3D, the linear term of 
q in Eq. ([8]) is now absent in 3D, and hence we could 
not exclude the possibility of Q being zero. In fact, the 
one-loop calculation at T = for a system in continuum 
shows that 



A 3D (g) = 




where x = q/2kp and d c — y 67r/(fci?m). As can be seen 
from Fig. |3(b)[ we found Q = within the one-loop 
approximation, implying a phenomenon of phase separa- 
tion for d > d c . As discussed above, this is sensitive to 
the short-range details of the interaction and no universal 
conclusion is available for case. 

Now we briefly discuss the multi-component case in 2D 
(e.g. the degrees of freedom of hyperfine spins) whose 
Hamiltonian is 

k.a 

+ EE^Ai^/' (14) 

k,k l ,q< 7 ' (J ' 



The stripe phase has a density modulation in the direc- 
tion perpendicular to the stripes, which can be detected 
directly via the measurement of the local densities. In ad- 
dition, any scattering experiments would show two inter- 
ference Bragg peaks at momentum ±Q in a stripe phase 
due to this density modulation. 

There are also other indirect tools to detect a stripe 
phase. Considering a density wave with density varying 
along the y axis , whose density profile (to the leading 
order) is 

P(x,y) = Po + PiCOs(Qy + (/)), (15) 

with (x, y) being the 2D coordinate in the real space, if 
we place this system in a shallow potential trap and then 
introduce an extra narrow ID potential well along the 
direction of the stripe, for example, 

V(x,y) = -aS(y- y ), (16) 

it is energetically favorable to have y = yo being a high 
density region due to the potential well V(x, y). We em- 
phasize that in reality, this ID potential well of V(x,y) 
shall have a finite width in the y direction. However, as 
long as the width is much smaller than the wavelength 
of the stripe, one can treat it as a 5-function to the lead- 
ing order approximation. If we move the potential trap 
along the x direction but keep the location of the po- 
tential well V(x,y) unchanged, the stripes will largely 
remain at their initial positions if the displacement of 
the trap is less than half the wave-length of the stripes 
Tt/Q, and the stripes will jump over a distance of 2tt/Q if 
the displacement becomes larger than ir/Q to same the 
potential energy. On the other hand, the jump would not 
happen in a homogeneous state. This jump will lead to a 
center of mass oscillation for the trapped particles due to 
the lack of dissipation, which could serve as an indirect 
signature for the stripe phase. 
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This technique is in close analogy to the "pinning" ef- 
fect, which was used to experimentally distinguish the 
nematic phase from a putative stripe phase in two- 
dimensional electron gases (2DEG) under high magnetic 
field [40|-|42|. There, the impurities in the crystal provide 
similar effects as the local potential in Eq. (|16[) . which 
will pin down the stripe if the stripe exists. In principle, 
the presence of stripes will leads to the pinning effect, 
which shall result in nonlinear signal in the I-V curve at 
low bias. However, in the 2DEG under high magnetic 
field, the absence of such a nonlinear I-V curve was in- 
terpreted as the absence of the stripe ordering. Different 
from the 2DEG, in the case we studied here, the simi- 
lar pinning effect would lead to oscillations for a system 
with stripe ordering, due to the absence of dissipation in 
ultracold dipolar gases. 

VI. DISCUSSIONS 

Since a real system has a finite temperature, the contri- 
bution of the thermal fluctuations comes into play. Their 
effects can be studied based on the symmetry-breaking 
pattern of the ordered phase. For a 2D system in the 
continuum, the homogeneous phase at d < dc has con- 
tinuous translational and two-fold rotational symmetries 
(rotation by tt) at 8 ^ 0. For d > dc, the stripe phase 
spontaneously breaks the continuous translational sym- 
metry in one direction and leads to one gapless Goldstone 
mode as required by symmetry. In a 2D system with 
d > dc, if we increase the temperature from T = 0, such 
a Goldstone mode will destroy the long-range order of 
the stripe phase. However, a quasi-long-range power-law 
correlation would still remain for temperature below a 
transition temperature Tkt [39]. For T < Tkt, although 
the real-space density oscillation of the stripes is hard 
to observe due to the lack of a true long-range order, 
the density profile in Fourier space will have two peaks 
located at ±Q which decays as a power-law function of 
l^f Q\ . Above Tkt, the correlation becomes short-ranged 
and the phase transition at T = Tkt belongs to the usual 
Kosterlitz— Thouless (KT) university class. 

For dipolar molecules in a 2D optical lattice, if the 
wavevector of the stripe phase is commensurate with the 
lattice wavevector, the stripe phase breaks no continu- 
ous symmetry and hence no gapless Goldstone mode is 
present. Therefore, the long-range order of the stripe 
phase remains at finite T below a critical temperature 
Tc- However, if the stripe phase is incommensurate with 
the underlying lattice, a KT transition is expected. 



In the study of Sees. Mil and IIV1 we were focusing on 
on the susceptibility shown in Eq. (JSJ) - Effectively, this 
procedure is equivalent and closely related to the liner- 
response approximation. Thus, the higher-order terms 
beyond the liner-response approximation were ignored in 
our study. These higher order terms are highly nontrivial 
and could become nonlocal at certain wave-vector. For 
example, the free energy contains terms proportional to 
p(q) 5 / 2 or p(q) 7 / 2 at q ~ 2k F , as shown in Refs. 
and [34J. Although these terms could be studied per- 
turbatively within certain approximations, there is no 
available technique to take care of them nonperturba- 
tively in a general way. In general, the contributions of 
these higher order terms are non-universal. In the case of 
a second-order phase transition, from the homogeneous 
phase to the inhomogeneous phase, these higher order 
term are unimportant at low energy and hence all the 
conclusions above remains. However, if the transition is 
first order, their contributions could become dominant. 
For a system with short-range interactions only, a first 
order transition is usually believed to leads to phase sep- 
aration (as in the liquid/gas-solid transition). However, 
it was also shown that a stripe phase was favored in the 
presence of long-range interactions like Coulomb 0, EE] 
or dipolar interactions [3q |. 

In summary, we find that the 2D dipolar Fermi gas 
undergoes a stripe ordering at large dipolar interaction 
strength, i.e., the density wave instability at finite mo- 
mentum, instead of phase separation. For the single 
component case, this instability occurs regardless of the 
dipole orientation. This arises from the fact that the 
Fourier transform of the dipolar interaction increases lin- 
early with wavevector. For the multicomponent case, due 
to the short-range Hartree interaction, the stripe insta- 
bility only exists at 8 > 9q. 

Note added 

After the completion of first online version of the 
manuscript, we learned the work by Yamaguchi et al. 
[46| , in which the charge-density- wave instability of dipo- 
lar Fermi gas are studied using mean-field approximation. 

Acknowledgments 

K. S. and S. D. S. are supported by JQI-NSF-PFC 
and AFOSR-MURI. C. W. is supported by NSF under 
No. DMR-0804775, and Sloan Research Foundation. 



[1] S. Ospelkaus, A. Pe'er, K.-K. Ni, J. J. Zirbel, B. Neyen- 
huis, S. Kotochigova, P. S. Julienne, J. Ye, and D. S. Jin, 
Nature Physics 4, 622 (2008). 

[2] K.-K. Ni, S. Ospelkaus, M. H. G. de Miranda, A. Pe'er, 



B. Neyenhuis, J. J. Zirbel, S. Kotochigova, P. S. Julienne, 
D. S. Jin,, and J. Ye, Science 322, 231 (2008). 
[3] A. Griesmaier, J. Werner, S. Hensler, J. Stuhler, and T. 
Pfau, Phys. Rev. Lett. 94, 160401 (2005). 



6 



[4] R. Chicireanu, A. Pouderous, R. Barbe, B. Laburthe- 
Tolra, E. Marechal, L. Vernac, J.-C. Keller, and O. Gor- 
ceix, Phys. Rev. A 73, 053406 (2006). 

[5] J. J. McClelland and J. L. Hanssen, Phys. Rev. Lett. 96, 
143005 (2006). 

[6] M. Lu, S. H. Youn, and B. L. Lev, Phys. Rev. Lett. 104, 
063001 (2010). 

[7] T. Koch, T. Lahaye, J. Metz, B. Frohlich, A. Griesmaier, 

and T. Pfau, Nature Physics 4, 218 (2008). 
[8] T. Lahaye, C. Menotti, L. Santos, M. Lewenstein, and T. 

Pfau, Rep. Prog. Phys. 72, 126401 (2009). 
[9] C. Menotti, M. Lewenstein, T. Lahaye, and T. Pfau, 

larXiv:0711.3422l (2007). 
[10] T. Sogo, L. He, T. Miyakawa, S. Yi, H. Lu, and H. Pu, 

New J. Phys. 11, 055017 (2008). 
[11] T. Miyakawa, T. Sogo, and H. Pu, Phys. Rev. A 77, 

061603(R) (2008). 

[12] S. Ronen and J. Bohn, larXiv:0906.3753l (2009). 

[13] C.-K. Chan, C. Wu, W.-C. Lee, and S. Das Sarma, Phys. 

Rev. A 81, 023602 (2010). 
[14] B. M. Fregoso, K. Sun, E. Fradkin, and B. L. Lev, New 

J. Phys. 11, 103003 (2009). 
[15] C. Lin, E. Zhao, and W. V. Liu, Phys. Rev. B 81, 045115 

(2010). 

[16] J. P. Kestner and S. Das Sarma. larXiv:1001.4763l (2010). 
[17] M. A. Baranov, M. S. Mar'enko, V. S. Rychkov, and 

G. V. Shlyapnikov, Phys. Rev. A 66, 013606 (2002). 
[18] M. A. Baranov, L. Dobrek, and M. Lewenstein, Phys. 

Rev. Lett. 92, 250403 (2004). 
[19] M. A. Baranov, Phys. Rep. 464, 71 (2008). 
[20] L. You and M. Marinescu, Phys. Rev. A 60, 2324 (1999). 
[21] G. M. Bruun and E. Taylor, Phys. Rev. Lett. 101, 245301 

(2008). 

[22] C. Zhao, L. Jiang and X. Liu and W. M. Liu and X. Zou 

and H. Pu. larXiv:0910.4775l (2009). 
[23] K. V. Samokhin and M. S. Mar'enko, Phys. Rev. Lett. 

97, 197003 (2006). 
[24] T. Shi, J. N . Zhang, C. P. Sun, and S. Yi, 

larXiv:0910.405Tl 2009. 
[25] C. Wu and J. E. Hirsch, Phys. Rev. B 81, 020508(R) 

(2010). 

[26] B. M. Fregoso and E. Fradkin. larXiv: 100 1.4 1671 2010. 
[27] B. M. Fregoso and E. Fradkin, Phys. Rev. Lett. 103, 



205301 (2009). 
[28] A. W. Overhauser, Phys. Rev. 128, 1437 (1962). 
[29] S. A. Kivelson, E. Fradkin, and V. J. Emery, Nature 393, 

550 (1998). 

[30] S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, 

J. M. Tranquada, A. Kapitulnik, and C. Howald, Rev. 

Mod. Phys. 75, 1201 (2003) 

[31] J.M. Tranquada, |arXiv:cond-mat/0512lT5| (2005). 

[32] E. Fradkin, S. A. Kivelson, M. J. Lawler, J. P. Eisenstein, 

A. P. Mackenzie. larXiv:0910.4166l (2009). 
[33] S. A. Brazovskii, Soviet Physics JETP 41, 85 (1975). 
[34] K. Sun, B. M. Fregoso, M. J. Lawler, E. Fradkin, Phys. 

Rev. B 78, 085124 (2008). 
[35] D. Amit, Field theory, the renormalization group, and 

critical phenomena (World Scientific, Singapore, 1984). 
[36] G. D. Mahan, Many Particle Physics (Plenum, New 

York, 1990). 

[37] The degeneracy between Q and —Q is due to the space- 
inversion symmetry r — > —r. 

[38] B. Spivak and S. A. Kivelson, Phys. Rev. B 70, 155114 
(2004). 

[39] P. M. Chaikin and T. C. Lubensky, Principles of Con- 
densed Matter Physics (Cambridge University Press, 
Cambridge, UK, 1998). 

[40] K. B. Cooper, M. P. Lilly, J. P. Eisenstein, T. Jungwirth, 
L. N. Pfeiffer, and K. W. West, Solid State Commun. 
119, 89 (2001). 

[41] K. B. Cooper, M. P. Lilly, J. P. Eisenstein, L. N. Pfeiffer, 

and K. W. West, Phys. Rev. B 65, 241313(R) (2002). 
[42] K. B. Cooper, J. P. Eisenstein, L. N. Pfeiffer, and K. W. 

West, Phys. Rev. Lett. 90, 226803 (2003). 
[43] D. L. Maslov, A. V. Chubukov, and R. Saha, Phys. Rev. 

B 74, 220402 (R) (2006). 
[44] V. J. Emery and S. A. Kivelson, Physica C 209, 597 

(1993). 

[45] S. A. Kivelson and V. J. Emery, in it Strongly Corre- 
lated Electron Materials: The Los Alamos Symposium 
1993, edited by K. Bedell, Z. Wang, D. E. Meltzer, A. V. 
Balatsky, and E. Abrahams (Addison- Wesley, Redwood 
City, 1994), pp. 619650. 

[46] Y. Yamaguchi, T. Sogo, T. Ito, T. Miyakawa, 
arXivl003.2029, (2010). 



